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ABSTRACT 

We have examined the interstellar scintillations of the pulsars in the double pulsar binary system. Near 
the time of the eclipse of pulsar A by the magnetosphere of B, the scintillations from both pulsars 
should be highly correlated because the radiation is passing through the same interstellar plasma. We 
report confirmation of this effect using 820 and 1400 MHz observations made with the Green Bank 
Telescope. The correlation allows us to constrain the projected relative position of the two pulsars 
at closest approach to be 4000 ± 2000 km, corresponding to an inclination that is only 0.29 ± 0.14° 
away from 90°. It also produces a two-dimensional map of the spatial correlation of the interstellar 
scintillation. This shows that the interstellar medium in the direction of the pulsars is significantly 
anisotropic. When this anisotropy is included in the orbital fitting, the transverse velocity of the 
center of mass is reduced from the previously published value of 141 ±8.5 km s""'^ to 66 ±15 km s^^. 
Subject headings: pulsars: general - pulsars: individual (J0737— 3039) - ISM: general - binaries: 
general 



1. INTRODUCTION 

The two pulsars in the recently discovered double 
pulsar binary system, "A" and "B", have periods of 
23 ms and 2.8 s, respectively. They are in a highly 
relativistic orbit with period of only 2.4 hrs and ec- 
centricity of 0.088 and have an average separat ion of 
only 2.9 It-s (|Burgav et al. 2003tlRTOe et al. 200|) . The 
phenomenology exhibited by this system is extremely 
rich. Pulsar A is eclips e d when it passes behind pul- 
sar B (|Lvne et al. 20041 iKaspi et al. 2004|) . Modula- 
tion of the eclipse with the period of B is consistent 
with the eclipse being caused by synchrotron absorp- 
tio n by the plasma surroun ding the magnetosphere of 
B (jMcLauehlin et al. 2004b|) . Throughout the orbit, 
the flux density of the B pulsar varies dramatically 
and systematically, indica ting significant interaction 
between the two pulsa rs l)Ramachandran et al. 200^ 
iMcLaughlin et al. 2004a|) . Only ~ 0.5 kpc away, and 
with orbital velocities reaching over 300 km s"\ this 
system is an ideal laboratory for studies of general rela- 
tivity, binary kinematics and, in the case of this paper, 
the binary orbital geometry and the interstellar medium. 

Pulse timing measurements provide many of the pa- 
rameters of the system, including the exact masses of 
the two pulsars, with exquisite precision. From measure- 
ments of the Shapiro delay, the incli nation of the sys - 
tem can be constrained to be 87° ±3° (|Lvne et al. 200^ . 
Measurements of the time scale of the interstellar scin- 
tillation (ISS) over an orbit can be used to estimate 
the orbital inclination, the longitude of periastron, and 
the transverse velocity of the center of mass of the sys- 
tem. This technique was first suggested by Lyne & 
Smith (198 2) and has been applied to PSR B0655±64 
()Lvne 1984|) and PSR J1141-6545 (Ord, Bailes & van 
Straten 2002). In general, there are two solutions for the 

^ University of California, San Diego 

^ Jodrell Bank Observatory, University of Manchester, Maccles- 
field, Cheshire, SKll 9DL, UK 

^ Massachusetts Institute of Technology, Haystack Observatory, 
Westford, MA 01886 



orbital inclination that fit the data equally well, but one 
may fit the expected mass of the pulsar better. In the 
case of the double pulsar system, the masses of A and 
B are already accurately determined through timing, re- 
moving this ambiguity. 

Ransom et al. (2004) (hereafter R04) applied the tech- 
nique to the double pulsar system, calculating an inclina- 
tion of 88.7° ±0.9° and transverse velocities of the center 
of mass parallel and perpendicular to the orbital line of 
nodes in the plane of the sky of Vpar — 96.0 ± 3.7 km 
and VJjorp = 103.1 ± 7.7 km s~^, respectively. This cal- 
culation relied on the implicit assumption that the tur- 
bulence in the interstellar medium (ISM) is isotropic. If 
the turbulence is anisotropic the transverse velocity es- 
timated by this method depends strongly on the axial 
ratio and orientation of the anisotropy and can be much 
lower than for isotropic turbulence. Note that neither 
this method nor timing can resolve the sense of the in- 
clination, i.e. whether the angular momentum vector 
points towards or away from the Earth. Furthermore 
neither method can resolve the direction of the parallel 
velocity with respect to the angular momentum vector. 

In this paper, we examine the correlations of the A and 
B pulsars' scintillations near the time of the A eclipse, 
i.e. when the lines of sight are closest to one another. 
We show how these correlations constrain the orbital ge- 
ometry more tightly than was previously possible and 
partially resolve the ambiguity in the sense of the incli- 
nation and the direction of the perpendicular velocity. 
We also show that the ISM is significantly anisotropic 
and that the transverse velocity of the center of mass is 
much smaller than was estimated with an isotropic model 
of the ISM. 

2. OBSERVATIONS AND ANALYSIS 

Observations of the double pulsar system with the 100- 
m Green Bank Telescope took place in December 2003 
and January 2004 under an "Exploratory Time" pro- 
posal as reported by R04. The 820-MHz data used for 
this analysis were acquired with the GBT Spectrometer 
SPIGOT card with a sampling time of 40.96 /xs on each of 
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1024 frequency channels covering a 50-MHz bandwidth. 
The 1400-MHz data were acquired using the Berkeley- 
Caltech Pulsar Machine (BCPM) using a sampling in- 
terval of 72 fj,s on each of 96 channels covering a 96-MHz 
bandwidth. A total of 5 hours of data at 820 MHz and 
6.1 hours of data at 1400 MHz were obtained. 

Our first step in the analysis of these data was to create 
dynamic spectra for each pulsar in each frequency band. 
In order to do this, we folded each frequency channel 
modulo the pulsar period using the ephemeris from Lyne 
et al. (2004) and freely available software tools (Lorimer 
et al. 2000). Individual frequency channels were shifted 
with respect to each other using a dispersion measure 
(DM) of 48.9 cm-3 pc (Lyne et al. 2004). Individual 
profiles were added so that the time resolution was 10 s 
at 1400 MHz and 5 s at 820 MHz. We formed on-pulse 
spectra by integrating over an "on-pulse" window and 
subtracted off-pulse spectra formed by integrating over 
an equivalent "off-pulse" window. We did not attempt to 
correct for variations in the system gain over the band- 
pass because such corrections increased the noise and the 
gain variations did not appear to bias the results. The 
dynamic spectra are essentially the same as those shown 
by R04. 

In order to estimate the time scale we first split the 
dynamic spectrum of A into segments of length 5.3 min 
and 10 min at 820 and 1400 MHz respectively, subtracted 
the mean from each spectrum, and then computed the 
temporal autocorrelation function p{t) of each segment. 
We then fit a theoretical model beginning at the sec- 
ond point to determine the time scale Atd, as defined 
in Cordes (1986). We used the theoretical form for a 
Kolmogorov scattering medium p{t) = exp[—{t/Atd)^^'^] 
rather than the commonly used Gaussian function p(t) = 
exp[— (t/A^d)^]- The rms error in the Atd estimates is 
much smaller at 820 MHz because the 820 MHz spectra 
have 10 times as many degrees of freedom as the 1400 
MHz spectra. 

Normally, scintillation is thought of as primarily a tem- 
poral process I{t), i.e. "twinkling" caused by the mo- 
tion of a spatial diffraction pattern /(r) past the ob- 
server with transverse velocity Vjss- If the diffraction 
pattern is isotropic with spatial scale 5*133, the time se- 
ries will have a time scale Tis3 = S'i3s/|Vi3s|. The 
spatial scale can be estimated from the bandwidth and 
the distribution of turb ulence along the line of sight 
(jCordes fc Rickett 1998j) . so that one may estimate the 
velocity as |Vis3| = 5133 /Tiss. The transverse velocity 
of the center of mass (Vpar, V^erp) and 5'iss are derived by 
fitting a model to the measured |Vi3s((/))|, where (j) is the 
orbital phase. As noted by R04 the result is calibrated 
relative to the precisely known orbital velocity of the A 
pulsar. However, if the ISM is anisotropic, 5*133 depends 
on the direction of V13S and the model must be modified 
accordingly (see SQJ- Our estimate of (V^ar,Vporp) us- 
ing an isotropic model matches that of R04 within their 
quoted errors. 

3. CROSS-CORRELATION OF PULSARS A AND B 

Near the time of the eclipse of pulsar A, the lines of 
sight to the two pulsars will pass through nearly the 
same ISM and their scintillations should therefore be 
correlated. The pulsars sample the ISM at their pro- 
jected locations rA{tA) = rA(0) + Va^a and rBits) = 



rB(0) -I- VetB which, of course, depend on the times of 
observations. Va and Vb are the vector sums of the 
transverse center of mass velocity (V^ar,V^erp) with the 
transverse orbital velocities Vao and Vbo- We can re- 
gard Va and Vb as essentially constant near the time 
of the eclipse. We define tA and ts with respect to the 
time of the superior conjunction of A (where A is at least 
partially eclipsed). Likewise, r\ and tb are defined with 
respect to the position of A at t^i 0. We define our 
coordinates as parallel and perpendicular to the line of 
nodes. Since the inclination is so close to 90 degrees, this 
is essentially parallel to pulsar B's orbital velocity. The 
location of B at superior conjunction is yB(0). The spa- 
tial cross correlation p{s = ta — tb) cannot be computed 
as a time average because the baseline s varies rapidly 
with time. However the averaging can be done over fre- 
quency because ISS also modulates the spectrum, and 
the spectral modulation is independent of velocity. 

In Figure n we present the cross-correlation p(tA,tB) 
of the dynamic spectra of A and B at 1400 MHz. This 
is an average of all three eclipses in the data set. The 
spectra were smoothed over three samples (30 s) in time, 
retaining 10 s sampling, in order to optimize the signal to 
noise ratio. The correlation is normalized by the square 
root of the product of the total variance from A and B, 
i.e. including noise. The marginal distributions show the 
fluxes of A and B as a function of time. The intensity of 
B is very low more than 200 s before the eclipse, mak- 
ing it impossible to correct the normalization for system 
noise over the range of (i^j^s) plotted in Figure ^ We 
corrected the correlation for noise in the region where the 
flux of B is significant (ts > —200 s) and fitted an ellipti- 
cal Kolmogorov model to the data. The peak correlation 
was 81% at Iap = 33.4±2.8 sec (i.e. after the eclipse) and 
tsp = — 22.7± 1.2 sec (i.e. before the eclipse). The noise 
correction is uncertain because the noise is not station- 
ary and we believe that the peak correlation is consistent 
with 100%. At this peak the A and B trajectories rA(iyi) 
and r-Bits) intersect as shown in Figure 2. 

The intersection equation rA(iAp) = ^B{tBp) yields 
Va^ap ~ VBi_BP + yB(0). The x component pro- 
vides a solution for the center of mass velocity com- 
ponent Vpar = (VAotAP + VBotBp)/{tAP - tsp)- Us- 
ing the known orbital speeds at the time of the eclipse 
{Vao = 329, Vbo = 352 km s~^ ), we directly obtain 
Vpar = 53 ± 17 km s^^ . This is clearly inconsistent with 
Vps^T- = 96 ± 4 km s~^ obtained by R04. We will show 
later that this inconsistency is due to anisotropic scatter- 
ing in the interstellar plasma. The y component of the 
intersection equation yields Vperp = yB(0)/(iAP — tsp)- 
As (<AP — ^Bp) is positive, we know that Vpcrp is in the 
same direction as yB(0), but we cannot determine if they 
are parallel or anti-parallel to the angular momentum. 
Using the value V^erp = 103 ± 7.7 km s~^ determined by 
R04 we find yB(0) — 5780 km. The trajectories rA(iA) 
and r-Bits) shown in Figure 2 are calculated using the ve- 
locities determined by R04 assuming isotropic scattering 
and this value of yB(0). The alternative track for r-Bits) 
shown without a dotted line corresponds to the ambigu- 
ity in the sense of the orbital inclination. The inclination 
is 0.41° from 90° if yB(0) is correct. However anisotropic 
scattering also modifies V^erp as we will show later, and 
this reduces yB(0). While the pulsars will again sample 
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played in Figure |31 Here wc have displayed a more lim- 
ited range than in Figure ^ and we have corrected the 
normalization of the correlations for system noise over 
the entire plot. This reduces the bias due to system 
noise but it creates a diagonal stripe of spurious corre- 
lation running from (0, -30000 km) to (60000 km, 0) 
where the flux of pulsar B is very low (so the correction 
factor is very large). When yB(0) is chosen correctly, 
the (noise-corrected) spatial correlation should be sym- 
metrical around the origin. The value of yB(0) which 
best centers the correlation is 4000 km. Note that be- 
cause the spatial scale of the ISS is larger than that 
of the eclipse, the spatial correlation is not badly dis- 
torted by A's eclipse. This apparent correlation is clearly 
anisotropic with an axial ratio near 2. However, since 
the mapping used the center of mass velocity obtained 
from an isotropic scattering model, this result is not self- 
consistent. Consequently, in 21 repeat the orbital 
time scale fitting including anisotropy. 



Fig. 1. — The cross correlation p{tA, is) averaged over the three 
ecHpses measured. The flux densities of A and B are shown as 
marginal distributions, showing clearly the eclipse of A and the 
strong orbital variation of B. The correlation has not been corrected 
for system noise. 



similar paths through the ISM when B passes behind A, 
the emission from B is very weak at this crossing and we 
cannot repeat this analysis for that part of the orbit. 
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Fig. 3. — The apparent spatial correlation p(sx, Sy) at 1400 MHz, 
where (xA,yA) and (xBiye) a^re the parallel and perpendicular 
components of and rg, respectively. The mapping velocities 
are those determined from an isotropic scattering model. An ellipse 
with an axial ratio of 2 with the major axis at -50° is drawn over 
the correlation to guide the eye. Here yeCO) = 4000 km is used 
because it provides a more symmetric result than the 5780 km 
shown in Figure 1^ 
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Fig. 2. — The trajectories of the A and B pulsars referenced to 
the location of the A pulsar at the time of its eclipse. The sample 
spacing is 10 s. The duration of the eclipse is indicated by open 
circles on the A trajectory, otherwise denoted by crosses connected 
by a dotted line. The plus signs connected by a dotted line rep- 
resent the B trajectory for an inclination of 90.41°, corresponding 
to an ye (0) of 5780 km. The plus signs with no dotted line show 
the B trajectory for the opposite sense of inclination. To illustrate 
the mapping more clearly, the baseline s = r 4(20) — ra(lO) where 
(t^jta) = (20, 10) is indicated by the double-ended arrow. 

The temporal correlation in Figure ^ can be mapped 
into a spatial correlation p{s = rj^{ty\) — r-B{tB)) if the 
velocities arc known. This will allow us to examine the 
isotropy assumption directly. The spatial correlation re- 
sulting from the trajectories shown in Figure |2] is dis- 



In Figure^ we present the temporal correlation at 820 
MHz, averaged over three eclipses. As at 1400 MHz, the 
spectra were smoothed over three sample intervals but 
the original 5 s sampling was retained. The spatial scale 
of ISS scales as A'^/" (eqn 44 of Rickctt, 1977), where 
a = 5/3 for a Kolmogorov scattering medium. Thus 
the scale at 820 MHz is smaller than at 1400 MHz by a 
factor of (820/1400)^'^ = 0.53, making it more compa- 
rable with the size of the eclipse. The raw correlations 
are lower than at 1400 MHz, but the estimation error is 
also lower. When the normalization is corrected for sys- 
tem noise the peak correlation reaches 100%. The spatial 
correlation at 820 MHz is shown in Figure |S1 using the 
same (Vpar, V^porp) as for Figure O The correlation is also 
anisotropic, but the axial ratio is somewhat smaller than 
at 1400 MHz. We attribute the differences to the greater 
distortion by the eclipse at 820 MHz. To make this cor- 
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relation symmetric about the origin we reduced yB(0) 
to 3000 km. It is possible that the difference between 
this value of yB(0) and that of 4000 km at 1400 MHz 
is due to differential refraction in the magnetosphere of 
B, but the difference is comparable with our estimate 
of the error. Allowing for this possibility we estimate 
that ye (0) = 4000 ± 2000 km which corresponds to an 
inchnation of 0.29° ± 0.14° away from 90°. 
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Fig. 4. — The average correlation pit Attn) for the three echpses 
at 820 MHz. The sample interval is 5 s. The origin is offset be- 
cause the first eclipse was very close to the start of the observation. 
The spatial scale of the ISS is much smaller than at 1400 MHz 
and is comparable to the scale of the eclipse, so that the correla- 
tion is heavily distorcd by the eclipse, and perhaps by propagation 
through the B magnetosphere. However, the peak remains in the 
lower left, as at 1400 MHz (Figure 1). 
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Fig. .5. — The spatial correlation p{sx, Sy) at 820 MHz. The edge 
of the transformation window can be seen because the first eclipse 
was very close to the start of the observations. This correlation 
has been corrected for 95% of the system noise. Full correction 
produced high errors during the eclipse. The ellipse drawn over 
the 1400 MHz data has been scaled down by a factor of 53% and 
over-plotted here. 



4. ANISOTROPY ANALYSIS FOR J0737-3039 



While the previous analysis showed that the ISS is 
anisotropic, the spatial mapping was carried out with ve- 
locities determined using an isotropic scattering model, 
hence the result is not self-consistent. In this section, we 
therefore modify the orbital fitting to include anisotropy 
and find a self-consistent spatial mapping. 

If the spatial correlation has elliptical symmetry, it can 
be written as p{s) = f{Q{s)), where Q{s) = as^ + bsy + 
cSxSy is the quadratic form which describes the ellipse. 
The form of / is less important; for example, for a Kol- 
mogorov power law, p(s) = exp[— ((5(s)/S'igg)^^^]. Here, 
S'iss is the 1/e spatial scale. The temporal correlation 
p{t) is determined by substitution of s = Vt, giving 
pit) = f{Q(V)t^), where V is the transverse velocity 
of the pulsar. The width Tigg of the temporal correla- 
tion is given by (9(V)Tjgg = Sfgg. In previous work, S'iss 
has been calculated from a theoretical scattering model, 
with the theoretical velocity Va(?!>) fit to the measured 
|Viss(?!')| = Siss/Tiss{4>) with an arbitrary scale factor 
(e.g. Ord et al. 2002). However, the preliminary calcula- 
tion of S'iss is not necessary, so we have simply included 
Siss in the fit instead of using a scale factor. The model 
we fit to the measured values of (1/Tiss) for pulsar A 
is [Q(Va)]°'^/5'iss, and the fitting parameters are Vpar, 
Vpcrp, and S'iss- The input parameters Ar and 9, where 
At is the axial ratio of the fitted ellipse and 9 is the 
position angle of the ellipse, must be given. The coeffi- 
cients of the quadratic form are: a = cos^O/ Ar + ArSin^ 9; 
b = 5111^9/ Ar + Arcos'^9; and c = 2sin6'cos6'(l/Ar - Ar). 
With this definition, Siss is the geometric mean of the 
major and minor axis scales. 

To find a self-consistent solution we use the 1400 MHz 
correlations, which are least distorted by the eclipse, with 
the 820 MHz measurements of 2iss(0) which have the 
higher signal to noise ratio. We then assumed an ax- 
ial ratio and orientation; computed the resulting center 
of mass velocity by fitting to the time scale vs orbital 
phase; used that velocity to create a spatial correlation; 
and estimated the axial ratio and orientation of that spa- 
tial correlation. When the correct axial ratio has been 
assumed it will be the same as that of the spatial correla- 
tion. We did this for a grid of axial ratio and orientation 
up to axial ratios of 10. To define a goodness of fit we 
used the technique suggested by Grail et al. (1997). We 
plotted the assumed anisotropy and that of the resulting 
spatial correlation on a modified Poincare sphere and 
used the distance between them as a measure of the er- 
ror. The best solution, near Ar = A and 9 = —15° with 
yB(0) = 4000 km, is shown in Figure 6. We find accept- 
able solutions for A,. > 3 with -20° < 9 < -10°. We 
searched Ar as large as 10, but we find it hard to believe 
that such axial ratios are feasible without corroborating 
evidence. The center of mass velocity is minimum at 
Ar = 4. If we restrict the axial ratio to 3 < < 6 
then 51 < |(V^ar, ^perp)| < 81 km s~^. The estimated 
Vpni- = 53± 17 km obtained earlier is consistent with 
this result and it may be possible with more data to use 
Vpar to improve the final model, but the data and analy- 
sis do not justify that at present. Thus we conclude that 
the center of mass velocity is much lower than previously 
thought. 

The best fit to the Tiss{(t>) measurements with this 
anisotropy has Siss = 9900 ± 300 km at 820 MHz. Note 
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that this is the geometric mean of the scales on the major 
and minor axes. This Siss corresponds to 18700 km at 
1400 MHz which agrees well with the spatial correlation 
shown in Figure 6. It is not possible to compare Siss with 
a theoretical calculation at present as the appropriate 
calculations for anisotropic scattering have not yet been 
done. 




XA — XB (10"' km) 



Fig. 6. — The final self-consistent spatial correlation p{sx,Sy) 
of the ISM in the direction of the J0737-3039 system for the ob- 
servations at 1400 MHz. The best-fit ellipse, with Ar = 4 and 
= —15°, is shown. 

Previous detections of anisotropic interstellar scatter- 
ing have come largely from angular broadening measure- 
ments of highly scattered objects. Measured axial ratios 
are in the range 1.3 to 3.0 (Wilkinson et al. 1994; Molnar 
et al. 1995; Desai & Fey 2001; Frail et al. 1994; Yuscf- 
Zadeh et al. 1994; Desai et al. 1994). An axial ratio 
of at least 4:1 was inferred from scintillation analysis of 
the quasar B0405-385 by Rickett at al. (2002). Though 
the latter result was identified as due to a discrete region 
only 15 - 30 pc from the Earth, it confirms that localized 
regions of the interstellar medium can cause axial ratios 
as high as we find here. 

It may be possible to confirm the anisotropy in this 
system with a more detailed analysis of the ISS using 
the "parabolic arc" phenomenon (Stinebring et al. 2001) 
and by extending the observations over a year so we can 
observe the annual modulation caused by the Earth's ve- 
locity. Parabolic arcs have been observed in this system 
and they provide a different view of the ISS. It is not 
fully independent of the time scale analysis but it has a 
different dependence on anisotropy which is the critical 
feature of our analysis. 

5. ANISOTROPY ANALYSIS IN GENERAL 

In fitting the Vigs('/>) measurements with anisotropic 
models, we realized that, because the inclination of the 
orbit is so close to 90°, Ar and 9 are completely de- 
generate parameters. The Viss(0) measurements them- 
selves can be fit equally well with any anisotropy (al- 
though with different Siss, ^par and Iporp)- The effect 
of anisotropy on Vp^r and V^orp is shown for Ar = i 
in Figure [7| The cause of the degeneracy can be de- 



duced from the equations for Va(0) directly (equations 
5-7 of Ord et al., 2002). The model for the degener- 
ate orbit as a function of orbital phase (f> has the form 
Va(0)^ = Acos20 -|- _Bsin0 + C (ignoring terms of the 
order of the eccentricity). Therefore, one can fit only 
three independent variables. Note that, in the case of the 
double pulsar system, the spatial correlations provide an 
estimate of the anisotropy and remove this degeneracy. 
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Fig. 7. — Vpar and Vperp derived from fitting of an anisotropic 
model to measurements of Vigg for Ar = 4 at various angles for 
J0737-3039. All fits had the same qu ality. The 95% confidence 
limits are shown. 

In general the Va(0)^ model has the form \^A{(f')'^ = 
Acos2(l) + Bsm2(j) + Csin0 4- Dcoscf) + E (again, ignoring 
terms of the order of the eccentricity). Thus one can fit 
five independent variables. In the system J1141— 6545, 
Ord et al. (2002) fit S'isS: ^par, V^erp, inclination, and 
the longitude of periastron. They found a most likely 
center of mass speed of 115 km s~^. We tested the effect 
of anisotropy on this system by fitting the measurements 
of Ord et al. (2002) and reproducing their result. Then, 
holding the longitude of periastron constant, we intro- 
duced an axial ratio and fit for the remaining parameters. 
We found a wide range of fits for 1.5 < A,. < 3, all of 
which were indistinguishable from the best isotropic fit. 
All of the fits we tried had greatly reduced Vpar and Vpcip. 
For example, with A^ = 2.1, we obtain = 0, inclination 
= 75°, Vpar = 20.5 km s"^ and Vpcrp = 17.1 km s"^ 

Although we do not know the value of the anisotropy 
for J1141— 6545, this is an important result, because it 
greatly weakens the case for high transverse velocities in 
binary systems. Observationally it makes the problem 
of measuring anisotropy in the ISM of more immediate 
interest. 

6. CONCLUSIONS 

The cross correlation between the ISS of the two pul- 
sars ill the double pulsar system has been measured. It 
tightly constrains the geometry of the eclipse. We es- 
timate the projected relative distance between the two 
pulsars at eclipse to be 4000 ± 2000 km. McLaughlin et 
al. (2004b) have shown that the eclipse duration varies, 
depending strongly on B pulse phase. However, if we take 
the lateral extent of the eclipse to be roughly 18,600 km 
(derived from the maximum 680 km s~^ relative veloc- 
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ity of the two pulsars and the average echpsc duration of 
27 s), the echpse must occur at an average radial distance 
from B of about 9000 km, only 7% of the 130,000 km ra- 
dius of the light cylinder of B. This is further evidence 
that the relativistic wind from A has blown away much 
of the magnetosphere of B. 

Our inclination estimate, 0.29° ±0.14° away from 90°, 
is consistent, within the errors, with the inclination de- 
termined from the Shapiro delay. However, the errors on 
the timing-derived inclination arc decreasing significantly 
with continued observations. If the estimate of inclina- 
tion determined through the ISS measurements differs 
from that measured from timing, we may be able to use 
the measurements of the ISS correlation at different fre- 
quencies to estimate the refraction of pulsar A by the 
magnetosphere of B near the eclipse and determine the 
density gradient in the magnetosphere. Note that the 
our inclination estimate is sufficiently small that gravi- 
tational lensing may be important (Lai & Rafikov 2004) . 
This effect would bias both the ISS and Shapiro delay 
estimates of inclination. 

The correlations show that the ISS in the direction of 
the system is anisotropic. When this anisotropy is in- 
cluded in the orbital analysis, the transverse velocity of 
the center of mass is reduced from 141 ± 8.5 km s~^ to 
66 ± 15 km s~^. We also see that a modest anisotropy in 
other binary systems, such as J1141— 6545, can greatly 
reduce the implied center-of-mass velocity. This de- 
creases the need to invoke large kick velocities and/or 
high precollapse core masses in double neutron star for- 
mation scenarios, as suggested by R04. In fact, as dis- 
cussed by Piran & Shaviv (2004), the low eccentricity of 
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the system, coupled with its location close to the Galac- 
tic plane, suggest a B progenitor mass less than 2 M©, 
and most likely around 1.45 Mq. This implies a non- 
standard, possibly white dwarf, progenitor. Willems & 
Kalogera (2004), who considered the orbital evolution of 
the system due to gravitational radiation and the orbital 
dynamics of asymmetric supernova explosions, likewise 
have difficulty explaining kick velocities of less than 60 
km s~^ given mass ranges derived by assuming a He- 
lium star progenitor and current models of Helium star 
evolution. Another way we may be able to confirm the 
low center of mass velocity of the system is through a 
measurement of the proper motion. Given the closeness 
of this system, such a measurement should be possible 
within a year. This will allow us to place better con- 
straints on the kick velocity, the progenitor mass of B 
and evolutionary scenarios, and will also allow us to test 
our main conclusion of a significantly anisotropic ISIVI. 



We are grateful to NRAO for making these data avail- 
able through the Rapid Response Science program and 
to Scott Ransom for taking the data. We thank Carlos 
Nava for his help estimating the time scales and fitting 
them with an anisotropic scattering model and Duncan 
Lorimer for writing the publically available SIGPROC 
package. The National Radio Astronomy Observatory 
is facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, 
Inc. Coles and Rickett acknowledge partial support from 
the NSF under grant AST 9988398. 



[McLaughlin et al. 2004b] McLaughhn, M. A. et al. 2004b, ApJ, 616, 
L131 

[Molbar et al. 1995]Molnar, L.A., Mutel, R.L., Reid, M.J. & 

Johnston, K.J. 1995, ApJ, 438, 708 
[Ord et al. 2002] Ord, S. M., Bailes, M. & van Straten, W. 2002, ApJ, 

574, L75 

[Piran & Shaviv 2004] Piran, T. & Shaviv, N. J. 2004, 

(astro-ph/0401553l 
[Ramachandran et al. 2004] Ramachandran, R. et al. 2004, ApJ, in 

press I astro-ph/0404392 I 
[Ransom et al. 2004]Ransom, S. M. et al. 2004, ApJ, 609, L71 
[Rickett 1977]Rickett, B. J., 1977, Ann. Rev. Astron. Astrophys., 

479 

[Rickett et al. 2002]Rickett, B. J., Kedziora-Chudczer, L., & 

Jauncey, D. L. 2002, ApJ, 581, 103 
[Stincbring et al. 2001] Stincbring, D. R. et al. 2001, ApJ, 549, L97 
[Wilkinson et al. 1994] Wilkinson, P.N., Narayan, R. & Spencer, 

R.E. 1994, MNRAS, 269, 67 
[Yusef-Zadeh et al. 1994]Yusef-Zadeh, F., Cotton, W., Wardle, M., 

Melia, F. & Roberts, D.A. 1994, ApJ, 434, L63 



